ON CONSTRAINED ANNEALED BOUNDS 
FOR PINNING AND WETTING MODELS 
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Abstract. The free energy of quenched disordered systems is bounded above by the 
free energy of the corresponding annealed system. This bound may be improved by ap- 
plying the annealing procedure, which is just Jensen inequality, after having modified the 
Hamiltonian in a way that the quenched expressions are left unchanged. This procedure 
is often viewed as a partial annealing or as a constrained annealing, in the sense that the 
term that is added may be interpreted as a Lagrange multiplier on the disorder variables. 

In this note we point out that, for a family of models, some of which have attracted 
much attention, the multipliers of the form of empirical averages of local functions cannot 
improve on the basic annealed bound from the viewpoint of characterizing the phase 
diagram. This class of multipliers is the one that is suitable for computations and it is 
often believed that in this class one can approximate arbitrarily well the quenched free 
energy. 



1. The framework and the main result 

1.1. The set— up (I): linear chain models. A number of disordered models of linear 
chains undergoing localization or pinning effects can be put into the following general 
framework. Let S := {Snj^^Q be a process with Sn taking values in Z*^, d G N := 
{1,2,...} and law P. 

The disorder in the system is given by a sequence oj := {uJn}^ of HD random variables 
of law P, with ujn taking values in F C M. As a matter of fact we could simply set F = M, 
however several examples that we will present deal with the case in which T is a finite set 
and in this situation our results require no measurability conditions. The disorder acts on 
the paths of S via an Hamiltonian that, for a system of size A^, is a function Hn^^^ of the 
trajectory S, but depending only on 5*0, Si, ... , Sn- One is interested in the properties of 
the probability measures Pn,lj defined by giving the density with respect to P: 

^^(S) = -^e^p{HN,^{S)), (1.1) 
dP Zn,w 

where Z^^^^ := E [exp {H^^^ {S))] is the normalization constant. Our attention focuses on 
the asymptotic behavior of log Z^^^. 

In the sequel we will assume: 
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Basic Hypothesis. There exists a sequence {Dn}^^ of subsets of Z such that P(5„, G 

TV — *oo 

Dn for n = 1,2, ... , N) x 1, namely 

hm ^AogP {SneDnioTn = 1,2,..., N) = 0, (1.2) 
and such that HN^i^{S) = if Sn £ Dn lor n = 1,2, . . . , N. 

One sees directly that this hypothesis implies 

lim inf 4 log Zn,u > lim ^ log P {Sn € D„ forn = 1, 2, . . . , A^) = 0, (1.3) 

P(dti;)-a.s.. We will assume that {(1/A^) log Zjv,tj}jY is a sequence of integrable random 
variables that converges in the (P(dti;)) sense and P(da;)-almost surely to a constant, 
the free energy, that we will call /. These assumptions are verified in the large majority of 
the interesting situations, for example whenever super /sub-additivity tools are applicable. 

Of course H1.3|l says that / > and one is lead to the natural question of whether / = 
or / > 0. In the instances that we are going to consider the free energy may be zero or 
positive according to some parameters from which Hi\f^;^{S) depends: / = and / > are 
associated to sharply different behaviors of the system. 

In order to establish upper bounds on / one may apply directly Jensen inequality 
{annealed bound) obtaining 



/ = hm i-E[logZTv,^] 

< lim inf 4 log IE[^^,^] =: 7g[0,oo] 



(1.4) 



and, in our context, if / = then / = 0. The annealed bound may be improved by adding 
to H]\f^i^{S) an integrable function : — > M such that E[AAr(a;)] = 0: in fact / as 
defined in the first line of (|1.4|) is unchanged by such transformation, while the second line 
of (|1.4() may depend on the choice of {An}]^. We stress that not only / is left unchanged 
by Hn,lu{S) —>■ Hiy^^{S) + Ap^{ui), but Pn,oj itself is left unchanged (for every N). Notice 
moreover that the optimal choice Aiy{uj) = — log Zjy^^ + E [log Zn uj] yields the equality in 

dn. 

In the sequel when we refer to / we mean that Zj^^^ is defined with respect to Hn^^^ 
satisfying the Basic Hypothesis (no An term added). 

1.2. The result. What we prove in this note is that 

Proposition 1.1. // / > then for every local bounded measurable function F : — > M 
such that E [F{uj)] = one has 

r / ^ 



lim inf — log EE 



n=0 



> 0, (1.5) 



where {9nUj)m = ^n+m- 



We can sum up this result by saying that when / = but /> it is of no use 
modifying the Hamiltonian by adding the empirical average of a (centered) local (bounded 
measurable) function. 
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Notice that requiring F{-) to be bounded and measurable is superfluous if F is a finite 
set. From now on the reader should read local as a short-cut for local, measurable and 
hounded. We take this occasion also to observe that in principle one should be able to 
extend the result in the direction of unbounded F{-) or of non IID disorder: this however 
requires additional assumptions and leads far from the spirit of this note. 

On a mathematical level it is not obvious that the free energy may be approximated 
via empirical averages of a local function of the disorder, because we are playing with 
an exchange of limits (recall the optimal choice of A^r above). But we remark that in 
the physical literature the approach of approximating the free energy via what can be 
viewed as a constrained annealed computation, the term "^^=0 ^i^n^) being interpreted 
as a Lagrange multiplier, is often considered as an effective way of approximating the 
quenched free energy. Here we mention in particular |2nj and jE] in which this point of 
view is taken up in a systematic way: the aim is to approach the quenched free energy by 
constrained annealing via local functions F that are more and more complex, the most 
natural example being linear combinations of correlations of higher and higher order. 

The proof of Proposition 11.11 is based on the simple observation that whenever A^r is 
centered 

^ log EE [exp {Hn,US) + An{oj))] > 

^logE[exp(A^H)] + ^logP(S„GZ)„forn=l,2,...,A^) =:Q^ + P^. (1.6) 

By hypothesis P/v = o(l) so one has to consider the asymptotic behavior of Qn. If 
liminfTvQAT > there is nothing to prove. So let us assume that liminfArQAr = 0: in 
this case the inferior limit of the left-hand side of (|1.6|) may be zero and we want to 
exclude this possibility when / > and A^iio) = J2n=o -^(.^n^)^ ^ local and centered (of 
course in this case limAr Qm does exist). And in Proposition 12. II below in fact we show that 
if logE [exp [A]\j(oj)y\ = o{N), then sup^^ |AAr(u;)| = o{N) and therefore the corresponding 
constrained annealing is just the standard annealing. 

Remark 1.2. We stress that our Basic Hypothesis is more general than it may look at 
first. As already observed, one has the freedom of adding to the Hamiltonian iJjv,Qj(5') 
any term that does not depend on S (but possibly does depend on to and N) without 
changing the model 'Pn,uj- It may therefore happen that the natural formulation of the 
Hamiltonian does not satisfy our Basic Hypothesis, but it does after a suitable additive 
correction. This happens for example in N1.2.3l below: the additive correction in that case 
is linear in uj and it corresponds to what in 21 is called first order Morita approximation. 
In these terms. Proposition II. II is saying that higher order Morita approximations cannot 
improve the bound on the critical curve found with the first order computation. 

Remark 1.3. In the Morita approach of |16| I2nj. when applied to spin systems, it was 
also taken for granted in the physics literature that the infinite volume measure describing 
the joint distribution of disorder variables and spin variables can be described as Gibbs 
measure with a proper (absolutely summable) Hamiltonian. This was shown to be false 
in general, and potentials with weaker summability properties are needed (3 EI- This 
phenomenon underlines from a different perspective that local dependence of the Morita 
potential on the disorder variables is not enough. 

Let us now look at applications of Proposition ll.il 
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1.2.1. Random rewards or penalties at the origin. Let S, Sq = G Z*^, be a random walk 
with centered IID non degenerate increments {Xn}n, {Xn)j £ {— li 0, 1} for j = 1,2, . . . , d, 
and 

N 

i^TV,.. = /3 J^(l + e^„)l{5„=o}- (1-7) 

n=l 

for /? > and £ > 0. The random variable uji is chosen such that E[exp(Awi)] < oo for 
every A € M, and centered. We write f{(3,e) for /: by super-additive arguments / exists 
and it is self-averaging (this observation is valid for all the models we consider and will 
not be repeated). We note that for e = the model can be solved, see e.g. |12| . and 
in particular 0) = if and only if /? < Pdd) ■= — log(l — P{S never comes back 
to 0)). Adding the disorder makes this model much more complex: the annealed bound 
yields /(/?, e) = if /? < (3c{d) — logE [exp(ea;i)] =: (3c- It is an open question whether 
Pc coincides with the quenched critical value or not, that is whether /(/?, e) = implies 
/3 < /?c or not. For references about this issue we refer to and see however also 
the next paragraph: the model we are considering can in fact be mapped to the wetting 
problem (0112). Proposition 11.11 applies to this context with Dn = {0}^ for every n jSl 
Ch. 3] and says that one cannot answer this question via constrained annealed bounds. 

1.2.2. Wetting models in 1 + d dimensions. Let S and lo be as in the previous example 
and 

rr _ jPEn=l (1 + ^^n) 1{(S„),=0} if (^n)^ > fom = 1, 2, . . . , iV 

I —oo otherwise. 

with /? > and e > 0. If one takes the directed walk viewpoint, that is if one considers 
the walk {(n, Sn)}n, then this is a model of a walk constrained above the (hyper-)plane 
Xd = and rewarded P, on the average, when touching this plane. If d = 1 then this is 
an effective model for a (l+l)-dimensional interface above a wall which mostly attracts 
it. As a matter of fact in this case there is essentially no loss of generality in considering 
d = 1, since localization is measured in terms of orthogonal displacements of the walk with 
respect to the wall and we may restrict ourselves to this coordinate. Once again if e = 
the model can be solved in detail, see e.g. ^21- Computing the critical P and deciding 
whether the annealed bound is sharp, at least for small e, is an unresolved and disputed 
question in the physical literature, see e.g. [Hlinilini- Proposition ll.ll applies with the choice 
Dn = Z'^-^ X N. 

1.2.3. Copolymer with adsorption models. For definiteness choose S to be a one dimen- 
sional simple random walk and take the directed walk viewpoint. Imagine that the space 
above the horizontal axis is filled with a solvent A, while below there is a solvent B. We 
choose LOi G {A, B} and for example 

N 

Hn%{S) = ^ (al{sign(5„)=+l,c^n=A} + ^l{sign(5n)=-l,L^„=B} + cl{5„=o}) (1-9) 
n=l 

with a, h and c real parameters and sign(S'„) = sign(5„_i) if 5„ = (this is just a trick 
to reward the bonds rather than the sites). In order to apply Proposition 11.11 one has to 
subtract a disorder dependent term, cf. Remark ll.2l if a > 6 we change the Hamiltonian 

N 

HnAS) ■■= H^iS) -Y,^hu.,.=A}- (1-10) 

n=l 
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without changing the measure Ptv.o; while the free energy has the trivial shift from / to 
/ — aP(ti;i = A). One can therefore choose Dn = Z"^"^ x N and Proposition 11.11 applies. 
This model has been considered for example in |21j . 

Note that if c = the model can be cast in a form that has been considered by a variety 
of authors (see e.g. [I3I21II1I1I1S1I171II3I3]): 

N 

HnAS) = a ^ K + /i) sign(S„), (1.11) 

n=l 

with uj taking values in M. Once again the Hamiltonian has to be corrected by subtracting 
the term A X]„(w„ + h) in order to apply Proposition ll.il One readily sees that H1.10() and 
(|1.11() are the same model when in the second case lo takes only the values ±1, A = +1 
and B = -1, and h = {a - b)/{a + b), A = (a + 6)/4. 

Proposition 11.11 acquires some interest in this context given the fact that the physical 
literature is rather split on the precise value of the critical curve and on whether the 
annealed bound is sharp or not, see 0| for details on this issue. In we present numerical 
evidence on the fact that the annealed curve does not coincide with the quenched one, and 
in view of Proposition II. II this would mean that constrained annealing via local functions 
cannot capture the phase diagram of the quenched system. 

1.2.4. Further linear chain models and observations. In spite of substantial numerical ev- 
idence that in several instances / = but / > 0, we are unaware of an interesting 
model for which this situation is rigorously known to happen. Consider however the case 
P(wi = +1) = P(a;i = -1) = 1/2 and 

N 

HnAS) = /? 5] (1 + eUrr) l{5„=n}, (1-12) 
n=l 

with P and e real numbers and S the standard simple symmetric random walk on Z. We 
observe that Proposition 11.11 applies to this case with Dn = {n}^ and that the model 
is solvable in detail. In particular /(/3, e) = {(3 — log 2) V 0, regardless of the value of e. 
The annealed computation instead yields /(/3, e) = {(3 + log cosh(e) — log 2) V 0. Notice in 
particular that the critical values of /?, respectively log 2 and log 2 — logcosh(e), differ as 
long as there is disorder in the system (e ^ 0). It is interesting to see in this toy model 
how the optimal choice of Af^f, mentioned at the end of S 11.11 is rather far from being the 
empirical average of a local function, when is large. 

Remark 1.4. We point out that we restricted our examples only to cases in which S is 
a simple random walk, but in principle our approach goes through for much more general 
models, like walks with correlated increments or self-interacting walks, see ^3 
example. And of course Sn takes values in Z'^ only for ease of exposition and can be 
easily generalized. Another important class of models to which our arguments apply is the 
disordered Poland-Scheraga one [lOj . 

1.3. The set— up (II): interface pinning models. It is natural to wonder whether one 
can go beyond the linear chain set-up. The answer is positive and we give the example of 
((i+ l)-dimensional effective interface models, d > 1, natural generalization of the (1 + 1)- 
dimensional interfaces considered in the previous section. By this we mean for example 
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the case of S := {S„}^^j^d with S'n G R and the law of 5 is P 



'N- 



P{d(f) OC exp j ^ U {(fn - fn')\ dfn 6o{d(pn) , (1-13) 



2 

n,n':\n—n'\=l I ngW 

where Vat = [-A^/2, n and [/(•) is a measurable function such limr^±oo U{r^ = 

+00 sufficiently rapidly to make the right-hand side of (|1.13|) integrable (note that we may 
assume U{-) to be even). As a matter of fact, in order to have a treatable model one has 
to restrict rather strongly the choice of U{-): interface models are extremely challenging 
even without introducing pinning potentials (or, of course, disorder). Connected to that 
is also the reason why we have chosen the continuous set-up for interface models: discrete 
models are even more challenging 13 . 

The disorder in the system this time is given by an IID field uj := ^-^d Hj^^^{S) 

depends only upon 5^ with n € Vn- ujq takes once again values in T. The definition (|1.1|) of 
P N,ui is unchanged and the Basic Hypothesis varies in the obvious way, that is we assume 
that there exists {^nlngz'* such that 

lim logP (S„ G Dn for n G IV) = 0, (1.14) 

and such that Hj\f^i^{S) = if G for every n G V/y. Like for linear chains we assume 
the existence of the quenched free energy, that is of the L^(P(da;)) and P(da;)-a.s. limit of 
the sequence |A~'^log Zj\[^^^j^ and like in the linear chain case we have < / < /, where 
/ is again the annealed free energy defined in analogy with ()1.4() . 

The punch-line of this section is that Proposition 11.11 holds in this new set-up and it is 
proven exactly in the same way: 

Proposition 1.5. /// > then for every local bounded measurable function F : F^'' — > R 
such that E [F{lj)] = one has 



liminf — 7 log EE 



exp Hn,US) + Yl ^(^^ 



ngAw 



> 0. (1.15) 



In order to give examples of applications we may consider the d + 1 dimensional model 
of random rewards and penalties near the origin, that is the case of 

Hn,u; = P Y1 (l+^^")l{5„e(-i,i)}, (1-16) 

neVjv 

but one can write natural straightforward generalizations of the wetting models and of 
the copolymer with adsorption. The Basic Hypothesis in all these cases is a probability 
estimate on what is known as an entropic repulsion event, that is, for example, the event 
that Sn > i for every n G Vn and one can for example show that such a probability is 
bounded below by exp (-ciV^'^) , c > 0, if [/(•) is and infrU"{r) > 0, see [ni and 
references therein. So in this case one may apply Proposition 11.11 to conclude that one 
cannot improve on the annealed bound by constraining via local functions. 

Two comments, of opposite spirit, are however in order (for details see the lecture notes 

my- 

(1) The Basic Hypothesis requires a substantially weaker estimate and it is reasonable 
to expect that one is able to verify it in greater generality. 
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(2) The understanding of the associated deterministic models (e = for random re- 
wards and wetting models and the annealed models in general) is still extremely 
partial. Somewhat satisfactory results are available for quadratic U{-), that is P is 
Gaussian, but even in this case one has to give up the precise estimates available 
for the linear chain case (like computing exactly Pc) and basic questions are still 
open. So the application of Proposition 11.51 while being relevant on a conceptual 
level, yields a result that has little quantitative content. 

2. On zero free energy and null potentials 

In this Section d > 1. Let {w„}„g2d be an IID family of random variables under the 
probability measure P, taking values in F = M. The law of uJi is denoted by z^. 

We are interested in the family A = {Ai\f}]\f^fq of empirical averages of a local function 
F, that is 

An{uj) = ^ F{9nuj) , (2.1) 

nSVjv 

where F : F^'' — > M depends only on the variables indexed by a finite set A C Z"^, that 
is F{uj) = F{uj') if ujn = oj'n for every n € A. Notice that, by standard (super-additivity) 
arguments, the limit 

L{F) := Ihn^ ^ logE [exp {A^iu;))] , (2.2) 

exists. Moreover, by Jensen's inequality, L{F) > E [F(a;)]. 
We will prove the following: 

Proposition 2.1. Assume thatE[F{uj)] = 0. If L{F) = 0, then 

lim supUAr(a;)| = 0. (2.3) 

Of course, since the result is uniform in uj, the proposition covers also the linear chain 
set-up, where one considers ^[Ar/2j+i^Af rather than Va^- 

Proof. We consider the potential, in the sense of ^3 Def. (2.2)], $ := {^^sj^^^d defined 

by 

\ F iO^n^) if there exists n such that OnB = A, 
I otherwise. 
Let V be the single spin reference measure jlU Def. (2.9)] and let us set 

Z%{u;) := J exp(/7*(a)) J] ^(^^n) 11 '^-.(^^n)' (2-5) 

with H^{a) := l^^.^nViv^e ^b{<^)- Note that An{-) differs from H^{-) only by boundary 
terms so that sup^ \A]\f{uj) — H^{uj)\ < CN'^~^ for some C > (we recall that F{-) is 
bounded). Therefore it suffices to show that (|2.,'-{j) holds with A]\[{-) replaced by H^{-). 

Let us consider the ^-invariant Gibbs measure // associated to the potential $, the 
existence of which is established in a standard way by taking infinite volume limits with 
periodic boundary conditions (if has unbounded support tightness follows from the 
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fact that -P(-) is bounded). By jlll Theorem (15.30)] the relative entropy density of 1^°° 
{iy°°{dLo) := rinGZ'* ^(d<^n)) with respect to fi exists and can be written as 

li'^oo W {''^ I ^) = W - / ^^-'^ 

where TYy^^ (i^°°|^) is the relative entropy of with respect to fi, when both measures 
are restricted to the a-algebra generated by the variables {uJn}neVN- We have of course 
used the standard definition of relative entropy, H{fii\iJ,2) = / log(d/ii/d/i2)d//i for //i 
and iJ,2 two probability measures with fii absolutely continuous with respect to //2- A last 
remark on formula (|2.6() is that it holds for any choice of iv: this is just the independence 
of the free energy on boundary conditions. This independence may be seen directly since 
log{Zff{u!) /Z^{u!')) = 0{N'^~^) uniformly in u! and oj' and this implies also that the first 
term in the right-hand side of (|2.6|) may be replaced by L{F). 

Notice now that both terms in the right-hand side of (|2.6|1 are zero, respectively by 
the hypotheses L{F) = and E[F(a;)] = 0, and therefore, as a consequence of the Gibbs 
variational principle jlll Theorem (15.37)], is a Gibbs measure with the same specifi- 
cation of /X, but of course is the Gibbs measure with potential ^^''^ identically equal to 
zero and single spin measure v. This means that $ — <I>('^)(= is a negligible potential, 
that is im Theorem (2.34)] the function 

(^b(^)-<^(^)) (2.7) 
does not depend on the variables u;„ for n € Vn- We can write 

B:BnVM^^ B-.BcVn B:Br\VN¥'%,B(lVN (2.8) 

=: /iv(u;) + iiAf('^), 

and since 7^* (w) does not depend on the a;„'s for n G Vn we may change in the right-hand 
side the configuration oj with oj defined by setting u;„ = a;„ for n G and a;„ = c, c an 
arbitrary fixed constant, for n E V^. Therefore, in random variable terms, we have 

rL%{uj) = CN + RNi^), (2.9) 

with cn = /at (5) (notice that it is not random and it depends only on the choice of c). 
From the immediate estimate sup^^ \R]\!{u))\ < CN'^~^ for some C = C{F) > it follows 
that for all lv 

CN-CN'^-^ < n%{uj) < CN + CN'^-\ (2.10) 
and the hypothesis L{F) = yields immediately lim7v-»oo cat/A^"^ = 0. Therefore 

sup|W^(a;)| < CN + CN'^-^ = o{N'^), (2.11) 
and the proof is complete. 
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